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The ability to wirelessly power electrical devices is becoming of greater urgency as a component of energy 
conservation and sustainability efforts. Due to health and safety concerns, most wireless power transfer 
(WPT) schemes utilize very low frequency, quasi-static, magnetic fields; power transfer occurs via 
magneto -inductive (MI) coupling between conducting loops serving as transmitter and receiver. At the 
"long range" regime - referring to distances larger than the diameter of the largest loop - WPT efficiency in 
free space falls off as (l/d) 6 ; power loss quickly approaches 100% and limits practical implementations of 
WPT to relatively tight distances between power source and device. A "superlens", however, can concentrate 
the magnetic near fields of a source. Here, we demonstrate the impact of a magnetic metamaterial (MM) 
superlens on long-range near-field WPT, quantitatively confirming in simulation and measurement at 13- 
16 MHz the conditions under which the superlens can enhance power transfer efficiency compared to the 
lens-less free-space system. 



The superlens, which can refocus not only propagating far-field waves but also non-propagating, near- 
field waves, has been one of the more provocative concepts to emerge from the field of metamaterials 1 . 
A superlens comprises a medium whose electric permittivity 8 and magnetic permeability |i both take 
on the value of —1. The superlens structure offers a means of controlling and manipulating the near-fields 
that would otherwise decay rapidly away from a source. Initially, the superlens was proposed in the context 
of optics, where its use was suggested as a means of forming an image with resolution greater than that 
implied by the diffraction limit 2 ' 3 . The superlens functions via the excitation of magnetic and electric surface 
modes that couple to the near fields of an object placed on one side of the slab, subsequently bringing them 
to a focus on the opposite side. Since electricity and magnetism are nearly decoupled in the near field, it was 
realized early on that a superlens with either 8 = — 1 or /n = —1 could focus the near field of electric or 
magnetic sources, respectively. Imaging with a superlens has been demonstrated at visible and infrared 
wavelengths using thin layers of materials such as silver or silicon carbide whose dielectric functions take 
the value of 8 = — 1 at particular wavelengths 4 . 

At low frequencies, where magnetism is much more prevalent, superlenses based on structured metamaterials 
characterized as artificial magnetic permeability media have been pursued for a variety of applications, including 
as flux guides to enhance resolution in magnetic resonance imaging 5,6 . More recently, as interest in WPT schemes 
has risen, the use of superlenses and other metamaterial-based components to enhance transfer efficiency has 
been suggested 7 . At the low frequencies typical of inductive WPT schemes, the excitation wavelength exceeds 
10 m, whereas the dimension of the coils and loops is on the decimeter or centimeter scale. Thus, the loops can be 
initially approximated as magnetic dipoles as a route to gaining an intuitive understanding of the limits and 
behavior of the WPT system. 

A general treatment of a magneto -inductive WPT system was carried out by Kurs et al., who made use of 
coupled-mode theory to calculate the expected efficiency of two self- resonant coils separated by a distance d in 
free space 8 . In that work and the subsequent analysis of Urzhumov et al. 9 , power transfer efficiency is defined as 
the power dissipated in a load placed on the receiver coil (Rx) divided by the total power dissipated in the 
transmitter (Tx) and receiver circuits, as well as in any intermediate "relay systems" (such as a repeater coil, 
lens or metamaterial layer): 
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Once the Rx and Tx circuits and an optional relay system are spe- 
cified, power dissipation rates can be calculated using the coupled 
mode theory for arbitrary Rx and Tx configurations 8 , or the simpli- 
fied coupled- dipole formalism for small resonant coils 9 . Coupled 
mode theory treats the Rx and Tx coils as resonators whose inter- 
action is indicated by a mode coupling coefficient, tc, whereas the 
coupled- dipole theory uses the conventional notion of mutual 
inductance, L 2 \. Both theories predict the same long-range behavior 
for the WPT efficiency (1), namely, rj ~ d~ 6 . The physical origin of 
this power law is due to the 1/d 3 dependence of magnetic field in the 
near field of a magnetic dipole source. The overall system efficiency 
in magneto-inductive schemes is therefore inherently limited by the 
divergence of magnetic flux in free space, because power transfer is 
ultimately related to the amount of flux from the first coil that can be 
captured in the "aperture" of the second coil. Given that any coil or 
loop will behave roughly as a magnetic dipole, it is unlikely that any 
redesign or engineering of the coils can possibly enhance WPT effi- 
ciency of any system that is already optimally impedance-matched. 

If the efficiency of a resonant WPT system is to be improved, a 
means must be found to recapture and refocus the otherwise diver- 
gent magnetic flux. Traditional lenses based on conventional materi- 
als only focus the far-fields, and therefore are irrelevant to this 
application. Near-field lenses of several configurations 10 " 12 have been 
proposed; however those configurations provide partial field focus- 
ing with efficiency significantly lower than the perfect focusing pro- 
mised by the superlens 13 . 

The use of a near-field superlens in conjunction with resonant 
power transfer was first considered by Wang et al 14 . Earlier, similar 
metamaterial lenses were proposed for magnetic resonance imaging 
applications 1516 . Examining both isotropic 14 as well as anisotropic 17 
versions of the negative-permeability medium, these works provided 
significant numerical evidence, and some experimental evidence, of 
enhanced efficiency. Here, we address the question raised by the 
studies in Ref. 14,17, as to whether metamaterial superlenses can 
improve WPT efficiency (and conversely, reduce power loss) in the 
long-range near-field transfer regime. 

Theoretically, a positive answer to this important question was 
given in Ref. 9. Figure 1 is based on the analytical solutions from 
Ref. 9, and it shows that a metamaterial with finite and realistic 
resistive/magnetic loss can help deliver magnetic flux and generate 
alternating current in a resonant Rx circuit, leading to overall 
improvement in WPT efficiency. This performance boost was pre- 
dicted in the long-range, high-load transfer regime, where both the 
transfer distance d and the resistive load in the receive circuit exceed 
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Figure 1 | Wireless power transfer efficiency of a D = 6 cm MM 
superlens (a layer with isotropic \i = —1 — ia) calculated using the 
coupled-dipole theory 9 . The solid curve is efficiency in free space. 
Resistances R Tx = 0.0 1Z 0 and Rrx = 0.1Z 0 . 




Figure 2 | Magnetic metamaterial and superlens assembly. (A) Double- 
periodic array of resonators generating magnetic response in x and y 
directions. (B) Overall view of a triple-layer superlens with the top z-facing 
sheet removed. (C) One of the small non-resonant coils used in experiments. 

certain thresholds. This Report presents experimental evidence that 
long-range WPT efficiency using a MM slab can exceed the max- 
imum efficiency obtainable, ceteris paribus, in free space. 

Results 

Negative-permeability superlens design. Geometrically, the super- 
lens is perhaps the simplest possible configuration for near-field 
focusing: it consists of a uniform layer of isotropic, negative perme- 
ability. This property simplifies the superlens design, which could be 
as simple as a double- or triple-periodic arrangement of identical unit 
cells. In addition, the superlens is translationally invariant, at least in 
the large aperture limit, which reduces the need for fine mechanical 
alignment of the WPT system relative to the lens. Figure 2 illustrates 
the geometry of our superlens implementation; the unit cell geome- 
try, design, homogenization, and effective permeability are detailed 
in Methods. The major challenge for designing negative-perme- 
ability metamaterials at relatively low frequencies (MHz regime 
and lower) is the design of sufficiently low loss, strong magnetic 
dipole resonators of deeply sub-wavelength dimensions. Here, we 
have opted for unit cells of size ~ 2 cm, and achieved sufficiently 
large inductance by using multi-turn planar coils 18 . 

Enhancement of Magneto-inductive coupling with finite aperture 
superlens. Once a metamaterial layer with desired effective perme- 
ability is designed, we can characterize its effect on the coupling 
between two magnetic dipoles. We approximate a finite-aperture 
slab as a disk of the same diameter, which enables highly efficient 
axisymmetric description along the lines of Ref. 9,17 in a 2D 
rotationally- symmetric geometry modeled in COMSOL Multiphy- 
sics. We use the retrieved permeability components of the meta- 
material slab, each fitted to a Lorentzian resonance shape, as the 
components of the slab diagonal permeability tensor /n = [fi T /n T 
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Figure 3 | Simulated (solid lines) and analytically calculated (dots) 
enhancement factor G = r p iab IT mc for transfer distance d = 16 cm and a 
slab with isotropic \i = [ujvJI/vUjv], anisotropic \i = [urHrM^v]> and 
indefinite |X = [1 1 \i N ] and \i = [HrHr 1]. Insets: Re(H 2 ) for the isotropic 
slab, revealing anti- symmetric and symmetric modes excited at the lower 
and higher-frequency enhancement peaks. 
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Figure 4 | Simulated (A) and measured (B) transmission coefficient on dB scale, 10log 10 (T), of a single-layer triple-negative metamaterial (MM), 
indefinite MM, and freespace for coil-to-coil distances of 8 cm (blue curves), 16 cm (green curves), and 24 cm (red curves). Simulated (C) and measured 
(D) lOlogxoCT) of a 3-layer MM and freespace for the same distances. 



u N ] (assuming all off-diagonal elements are zero and relative 
permittivity is e r = +1). The transmitting magnetic dipole is 
simulated as a current loop of radius r = 1 cm carrying a fixed 
current I = 1A. We define the magnetic field transmission 
coefficient as T = IH^/H^J 2 , where H Tx = I/2r is used to approxi- 
mate the total field radiating from the Tx loop, and H Rx is the field 
measured at a point on the axis -of- revolution a distance d away from 
the loop. It is convenient to define transmission enhancement factor 
due to the slab, using the ratio of transmission coefficients with and 
without the slab G = T lab IV ac . 

Since the simulated coil is small relative to the slab and can be 
approximated as a dipole, we use the simulation results to validate the 
analytical solutions for the inductance between two dipoles with and 
without the MM slab (L s jf and L™ c , respectively) as expressed in Ref. 
9. In Figure 3 we plot the simulated and analytical enhancement 

factors (the latter, defined as G= \U^ h / V^f , we evaluate numer- 
ically with Mathematica). Results indicate that the slab enhances the 
field transmission coefficient T by roughly a factor of five over an 
appreciable frequency region. In addition, we observe a striking 
agreement between simulation and analytical results as a function 
of frequency. 

Measurements of enhancement with non-resonant coils. To ex- 
perimentally verify the effect of the metamaterial "superlens" on 
WPT between two coils, we begin by constructing a slab from the 
fabricated MM. To assemble a 1 -layer slab, we first designed the x- 
and y-oriented MM elements in rows with opposite-facing slits, 
which can be assembled into a wine-crate pattern, as shown in 
Figure 2A. The z-oriented resonators were fabricated across a 
single sheet, which was then placed perpendicular to the wine- 
crate and held together by additional slits cut around the wine- 
crate. In a similar fashion, we assembled a 3-layer MM slab 
(Figure 2B). To demonstrate a dipole- to-dipole WPT system, we 
used two small copper coils with cross-sectional wire diameter of 
1.6 mm wound into a circular loop of radius 1 cm, aligned the 
coils coaxially and separated by distance d along their shared axis, 
and connected them to ports 1 and 2 of an Agilent Vector Network 
Analyzer (VNA). Since our coils are non-resonant, to compute the 
experimental WPT efficiency from the measured S -parameters we 
must first remove the losses resulting from the mismatch between the 
network analyzer and the coils. Calculations needed to account for 
these losses are detailed below in Methods. 

Figure 4 shows the simulated and measured transmission coeffi- 
cients for free space and in the presence of a 1- and 3-layer slab, 
across several coil-to-coil distances. In addition, for the 1 -layer slab, 
we report the transmission coefficient of the two indefinite MM 
layers: one composed of only the z-facing MMs while the other 
consists of the xy-intersecting MMs. 



Substrate 



_ ^OOum 




Figure 5 | (A) Schematic of the metamaterial unit cell suitable for Printed 
Circuit Board fabrication process. The coils on the opposite sides of the 
substrate are rotated with respect to each other and connected by vias. (B) 
Each resonator has 17 turns on each side of a 10 mil substrate. Metal line 
widths are 200 microns, and the gap between neighboring cells is 1 
millimeter. (C) The complete unit-cell consists of three mutually 
orthogonal resonators. (D) The 3-layer slab is composed of the unit-cells 
stacked along z. 
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Figure 6 | (A) Transverse component of the effective permeability fora 1 -layer MM slab, calculated using the conventional S-parameter retrieval methods 
(circles) and the field averaging method (solid line). With the latter method, we can also calculate the normal component (crosses). Blue and orange 
represent the real and imaginary parts of (a, respectively. (B) The transverse and normal component of a 3-layer slab, retrieved using the Field Averaging 
method. 



Discussion 

Our central experimental results as shown in Fig. 4 are generally in 
very good agreement with numerical simulations, with the exception 
of the curves corresponding to power transfer distance d = 8 cm. 
The latter case needs additional discussion, considering that the 
distance from Tx and Rx coils to the surface of the metamaterial slab 
(of thickness L = 6 cm) is only (d — L)/2 = 1 cm, or roughly one half 
of the lateral unit cell dimension (array periodicity). In our numerical 
simulations (Fig. 4A,C), we first simulate a single period of the 2D 
array, and retrieve all three principal values of the effective permeab- 
ility tensor. Then, we replace the slab with a layer of homogeneous 
magnetic medium having the same complex permeability compo- 
nents as retrieved. In doing so, near field effects not describable in 
terms of effective magnetic permeability are lost. The physical fields 
on the surface of an array can be Fourier transformed, and all waves 
with transverse wavenumbers k > k 0 = w/c are necessarily evan- 
escent. However, out of this infinite spectrum, the effective medium 

description is only adequate for components with k < ksioch = \ 



where a is the metamaterial lattice constant. The higher- order 
Fourier components with k > k B i 0C h decay rapidly away from the 
interface - as exp( — k B i och z) or faster. They can be picked up by a near 
field probe placed at z = l/k B i och (or closer), which corresponds to 
our measurements taken with d = 8 cm. However, these harmonics 
are negligibly small at distances d > 15 cm, which explains the 
excellent agreement between effective medium models and measure- 
ments in that regime (Fig. 4). 

To summarize, we have demonstrated that a resonant array acting 
as an effective medium with negative- definite magnetic permeability 
enhances near-field transmission of quasi-static magnetic fields 
between non-resonant magnetic loop antennas. Significant enhance- 
ment is seen only at frequencies where at least one component of 
effective magnetic permeability has a negative real part. The enhan- 
cements due to different components of permeability tensor are 



investigated (Fig. 3), and it is shown that the strongest effect is 
obtained when all three components of Re([i) are negative (see 
Fig. 3, red and blue curves). We attribute the latter effect to the 
excitation of magneto -inductive surface waves existing in triple- 
negative permeability layers 19 . Enhancements in power transmission 
coefficient in the range of + 15 to + 30 dB are observed for all transfer 
distances from 8-24 cm; those transfer distances are 4-12 times 
greater than the diameter of both transmitter and receiver coils. 
Much higher enhancement is anticipated with larger- aperture super- 
lenses and multi-turn, self-resonant Tx/Rx coils. 



Methods 

Magnetic metamaterial design and fabrication. Our aim is to design and fabricate 
an isotropic metamaterial (MM) exhibiting Re{/4 < 0 at 13.56 MHz with minimal 
losses. Assuming operational frequency close to 13.56 MHz, the free- space 
wavelength is X 0 ~ 22 m; a conventional metamaterial whose elements are roughly 
/l o /10 in size would be far too large for practical implementations. Instead, we must 
demonstrate the desired behavior with elements whose size is only several 
centimeters, on the order of /l o /1000. To achieve this, we significantly increase the 
metamaterial unit cell inductance by utilizing the double- sided rotated coil design 
shown in Figure 5, which sandwiches a substrate between two via-connected multi- 
turn coils. The coils are rotated with respect to one another to form a composite 
circuit in which the inductances of two individual coils are added in series, resulting in 
the total inductance improved by a factor of four relative to the configuration with 
inductance in parallel. 

We form the complete metamaterial unit cell by positioning three identical reso- 
nators perpendicular to each other, as shown in Figure 5C, and iteratively tweak the 
design in CST microwave studio using the standard S-parameters retrievals 20,21 to 
obtain Re{/4 = — 1 in the desired ISM frequency band. 

The figure of merit (FOM) for our design is the inverse loss-tangent ratio at the 
frequency where Re{fi\ = — 1: 



Re(fi) 



Im(fi) 



Re(n) = 



(2) 



To minimize the loss tangent we choose a low-loss Rogers 4350 substrate and 
construct the metamaterial using a 1 -ounce (34 um thick) copper clad; the skin depth 
in Cu at 10 MHz is about 20 um. The S-parameter retrieval method lets us compute 
the transverse permeability components \i x and \i y by setting periodic boundary 
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Figure 7 | Pre- measurement calibration moves the VNA reference planes to the end of the cables that are connected to the non-resonant coils. 



conditions along the x and)/ directions and enforcing a normally (z-) incident, 
transversely polarized plane wave. The final design, whose retrieved permeability is 
shown in Figure 6, utilizes coils with 17 turns on each side of the substrate. Each turn 
is 200 urn wide and the gap between turns is set to 200 um as well. To reduce the 
design's sensitivity to fabrication errors, we insert three vias 200 um in diameter into 
the outer- most leg in each coil, and increase the width of that leg to 500 um such that 
there are 150 um between each via and the metal's edge. The total unit-cell size, 
including a 1 mm gap between adjacent unit-cells, is 1.894 cm. 

Although all three orthogonal coils in the unit-cell are identical, this does not mean 
the MM's permeability is isotropic because the normal (z) and transverse (x,y) 
boundary conditions observed by fields propagating through the slab are significantly 
different from each other. 

Field-averaging homogenization method for finite-thickness, anisotropic 
metamaterial layers. While the standard S-parameter retrieval techniques 20 ' 21 allow 
one to compute the components of effective permeability and permittivity 
tangential to the surface of a MM layer, the normal components are difficult to 
retrieve since they are not excited by a normally incident, transversely polarized 
plane wave. Here we present a quasi- magneto static field averaging retrieval 
method and outline how we used it to compute both the transverse and normal 
components of our MM. Quasi- electro static field-averaging homogenization was 
described in detail in Ref. 22., where effective s was expressed through the 
capacitance of a unit cell submerged into a curl-free electric field. Here, we extend 
this method to quasistatic permeability retrieval using the electric- magnetic 
duality theorem. For brevity we assume the medium to be uniaxial with 
permeabilities fi T and fi N ; the method is applicable to a general orthotropic 
medium with an orthorhombic lattice. 

We simulate a unit cell of dimensions a x X a y X a z in COMSOL Multiphysics's RF 
module. Air surrounded the cell along z while periodic boundary conditions (BCs) are 
along x and y such that unit-cell is part of a slab. Across the faces normal to the z-axis 
we assign an Electric Field which varies with z, E x (z) =xE 0 exp(jk z z). By using 
/c z <Cl/a z , we can be sure that the electric field has a linear variation in z, and thus its 
curl is virtually uniform in the entire domain. From Faraday's law, this E-field leads to 
j>-polarized, uniform magnetic field H y . By using the duality theorem together with 
the definition of capacitance, C = sA/d = Q/V, we can replace s with ft, electric charge 
Q with magnetic charge Q m , and electric voltage V with magnetic potential V m , and 
extract effective permeability according to 



Qm 
Vm 



Qm d 
' V m A 



BdA 
Bdl 



(3) 



With H polarized along y, we obtain fi y , one of the permeability principal values, by 
substituting d = a y , A = a x X a z , and B = B y into (3). 

To compute fi x , the remaining transverse component, one replaces the incident 
field E x (z) with E y (z) and utilizes the appropriate fields and dimensions in (3). To 
compute the normal component fi z , however, an additional subtle change has to be 
made. We begin by exciting an incident field E y (x) of the form E y (x) =yE 0 exp (jk^x) , 
which gives rise to an almost uniform H-field H z . However, such a field still violates 
the periodic condition along x slightly; therefore on the x = const faces we use Floquet 
(phase- shifted periodic) boundary condition with the phase shift given by kxa x . We 
then compute \x z from (3), substituting d = a z , A = a x X a y . 

We perform field- averaging retrievals across the 8 MHz- 16 MHz frequency range 
and compute the transverse and normal components of ft for both the 1- and 3 -layers 
configurations. We then fit each retrieved permeability to a Lorentzian curve defined 
as 



ji=l- 



Fo) z 



Co 2 — (Dq +jcoy 



(4) 



where F is a constant representing the oscillator's strength, co 0 = 27i/is the angular 
resonance frequency, andy = ooJlQ. The resulting fitted parameters of the 1 -layer 
MM are 



Qli 

fo,N 



0.26 
= 18.5 
= 13 MHz 
= 0.41 



Fj = .39 
Cft = 17.7 
fl T = 11.8MHz 
o\ = 0.27 



(5) 



and the parameters fitted from the 3-Layer MM retrieval were calculated to be 

F^ = 0.35 F 3 T = .36 

Q 3 N = 17.5 Q 3 r = 16 

fl N = 12.25 MHz fl T = 12.24 MHz 

<7 3 =0.32 (7 3 r = 0.32 



(6) 



Here we have used the superscripts 1 and 3 to distinguish between the 1- and 3 -layer 
slabs, respectively, and the subscripts N and T to distinguish between the normal and 
transverse permeability components, respectively. Not surprisingly, comparing the 
fitted parameters for the 1 - and 3 -layer MM suggests that as more layers are added the 
metamaterial behaves in a more isotropic fashion. 

We note that the Lorentzian parameters provide a quality fit for the complex 
permeability curve only in the vicinity of the fundamental magnetic resonance 
studied. Good quality of fit is maintained at least through the frequency where 
Re(u) crosses zero (roughly 16 MHz), that is, in the entire frequency band of interest. 

Maximum transducer power gain calculations. Before we conduct WPT 
measurements with the non-resonant coils setup, we perform a calibration that moves 
reference planes of a VNA to the end of the cables that are connected to the coils (see 
Figure 7). This enables us to retrieve the direct coil-to-coil transmission efficiency. In 
Ref. 23. Pozar describes a suitable metric called the Maximum Transducer Power 
Gain, G Tmax , which is the gain that would be achieved if a lossless matching network 
was inserted between the NA's reference planes and the non-resonant loops. Pozar 
defines Gj max in terms of only S parameters; here we summarize the calculations 
outlined in Ref. 23. 

Transducer Power Gain G T is the ratio of power delivered to the load, P L , to the 
power available from the source, p s : 



p L _ |s 21 | 2 (i-|r s | 2 )(i-|r L | 2 
Ps |i-r sr j 2 |i-s 22 r i | 2 



(7) 



where T L = {Z L — Z 0 )/{Z L + Z 0 ) is the reflection coefficient seen looking toward the 
load, r s = (Z s — Z 0 )/(Z S + Z 0 ) is the reflection coefficient seen looking toward the 
source, and T in is the reflection coefficient seen looking toward the input of the two 
port network 



i-s 22 r L 



(8) 



where Z in is the impedance seen looking into port 1 of the terminated network. 
Similarly, Y out is the reflection coefficient seen looking into port 2 of the network 
when port 1 is terminated by Z s : 



i-s n r s 



(9) 



The Maximum Transducer Power Gain, G Tmax , occurs when T in = T* and T out = T* L . 
In the general case with a bilateral two port network, T in is affected by T out , and vice 
versa, so that the input and output must be matched simultaneously. Equating T* s and 
F* L with the RHS of (8) and (9), respectively, yields 



rs=s? 1 + 



(10) 



SCIENTIFIC REPORTS | 4:3642 | DOI: 1 0.1 038/srep03642 



5 



_S 22 -AT S 
' 1-SnTs 



(11) 



where A = S n S 22 ~ S 12 S 2 i. Substituting (10) into (9) and rearranging the terms 
results in the quadratic equation 

T s 2 (S n - AS* 22 ) + T s (| A| 2 - |S n | 2 + |S 22 1 2 - 1) + (S* n - A*S 22 ) = 0 (12) 

yielding the solutions 

r _ B 1 ±^Bj-4\C l \ 2 



2Q 



B 2 ±^B 2 2 -4\C 2 \ 2 
2C 2 



(13) 



(14) 



where ^ = 1 + |S n | 2 - |S 22 | 2 - |A| 2 , d =S n -AS* 22 ,B 2 = 1 + |S 22 | 2 - |S n | 2 - |A| 2 , 
and C 2 =S 22 —AS* U . To convert the recorded S-parameters to G Tmax , then, one must 
first calculate (in order) A, Q, C 2 , B 1? B 2 , r L , T s and r in and use these parameters in 
(12). 
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